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ABSTRACT 

A complete sample of 7.7 x 10^ galaxies with five-band imaging and spectroscopic redshifts from the 
Sloan Digital Sky Survey is used to determine the fraction of the optical luminosity density of the Local 
Universe (redshifts 0.02 < z < 0.22) emitted by red galaxies. The distribution in the space of rest- 
frame color, central surface brightness, and concentration is shown to be highly clustered and bimodal; 
galaxies fall primarily into one of two distinct classes. One class is red, concentrated and high in surface 
brightness; the other is bluer, less concentrated, and lower in central surface brightness. Elliptical and 
bulge-dominated galaxies preferentially belong to the red class. Even with a very restrictive definition 
of the red class that includes limits on color, surface brightness and concentration, the class comprises 
roughly one fifth of the number density of galaxies more luminous than 0.05 L* and produces two fifths of 
the total cosmic galaxy luminosity density at 0.7 fim. The natural interpretation is that a large fraction 
of the stellar mass density of the Local Universe is in very old stellar populations. 

Subject headings: cosmology: observations — galaxies: elliptical and lenticular, cD — galaxies: 
evolution — galaxies: fundamental parameters — galaxies: stellar content 



1. INTRODUCTION 

The local Universe sparkles with bright blue stellar pop- 
ulations in disk and irregular galaxies. However, both age 
and dust (circumstellar or interstellar) make a population 
dimmer and redder, so red populations contain much more 
stellar mass per unit luminosity. This makes it remarkable 
that the most luminous (visible and near-infrared) galaxies 
in the Universe are red galaxies (eg, Blanton et al 2001); 
their dominance in luminosity represents enormous domi- 
nance in stellar mass. 

The typical luminous red galaxies of the local Universe 
are kinematically hot (ie, not rotationally supported), have 
smooth radial profiles, have a narrow range of stellar mass- 
to-light ratios, show high a-element metallicities, are (rel- 
atively) iron-poor, and reside preferentially in the Uni- 
verse's most over-dense environments (as reviewed by Ko- 
rmendy & Djorgovski 1989, and Roberts & Haynes 1994). 
They are red because their stellar populations are old — 
comparable in age to the Universe itself — not because they 
are obscured by dust. 

It is now well established that the average cosmic star 
formation rate per comoving volume has been declining 
from redshifts near unity to the present day (eg, Pei & 



Fall 1995, Lilly et al 1996, and many more references cited 
in Hogg 2001). This observed trend is strong enough that 
the majority of the stars in the local Universe are expected 
to be old (> 6 Gyr). This appears consistent with observa- 
tions of the distribution of stellar ages at the present day 
(eg, Fukugita et al 1998). What is not currently agreed- 
upon is the detailed stellar age distribution at old ages. 
Direct measurements of the cosmic star formation rate at 
very early times are hindered by selection effects and the 
strong infiucnce of dust; measurements of age distributions 
of nearby stellar populations are confused by the asymp- 
totic similarity of spectral energy distributions with in- 
creasing age. 

In this Article we look at the luminosity density in red 
galaxies in the z ^ 0.1 local Universe using observations 
from the Sloan Digital Sky Survey (SDSS). We begin by 
suggesting a new method for selecting red galaxies de- 
signed to avoid contamination by young galaxies reddened 
by dust. We find that the red galaxies contribute a large 
fraction of the total galaxy luminosity density of the Uni- 
verse at A 0.7 fim, such a large fraction that they make 
up the majority of the cosmic stellar mass density. 

In what follows, a cosmological world model with 
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(f2M)^A) = (0.3,0.7) is adopted. The model has total 
age 9Ah~^ Gyr, where h is the Hubble constant in units 
of 100 kms^i Mpc-^ h = 0.67 makes the total age 14 Gyr 
(eg, Hogg 1999). 

2. THE SAMPLE 

2.1. Imaging and photometry 

The SDSS is taking u, g, r, i and ^;-band drift-scan 
imaging of the Northern Galactic Cap, and taking spec- 
tra of roughly 9 x 10^ galaxies in that region, most with 
r < 17.7 mag in the AB system (York et al 2000). The 
photometric system and imaging hardware are described 
in detail elsewhere (Fukugita et al 1996; Gunn et al 
1998; Smith, Tucker et al 2002). An automated image- 
processing system detects astronomical sources and mea- 
sures many photometric properties (Lupton et al 2001; 
Stoughton et al 2002, Pier et al 2002, Lupton et al in 
preparation). The photometric parameters of interest for 
this study are as follows: the Petrosian (1976) radius 6'pct 
is the angular radius at which the mean surface brightness 
of the source in the SDSS r-band image inside that radius 
is five times higher than the mean surface brightness in 
a narrow annulus centered on that radius. The Petrosian 
fluxes (with corresponding magnitudes mpd) are the total 
fluxes in the five SDSS bandpass images within a circu- 
lar aperture of radius 2 9r,pet (twice the Petrosian radius 
in the r-band). The half-light and 90-percent radii 
and 9go are the angular radii within which 50 and 90 per- 
cent of the Petrosian flux is found. All magnitudes are 
extiction-corrected with IRAS-based dust maps (Schlegel 
Finkbeiner & Davis 1998). 

The median magnitude of a galaxy in our sample is 
r = 17 mag, so it is detected in the imaging data at a 
signal-to-noise ratio S/N > 300; the photometric proper- 
ties of the galaxies in the sample are measured with high 
precision. The accuracy is limited by systematic errors at 
the few percent to ten percent level (Stoughton et al 2002). 

2.2. Spectroscopy 

SDSS "Main Sample" galaxies are chosen as spectro- 
scopic targets on the basis of r-band magnitude, angu- 
lar size (to remove stars), half-light surface brightness (to 
ensure spectroscopic success), and some positional con- 
straints from the finite size of the fibers in the fiber-fed 
multi-object spectrographs (Strauss et al 2002, Blanton et 
al 2002a). 

The SDSS spectra are at resolution 2000 with wave- 
length coverage 3900 < A < 9000 A. Typical spectra 
have signal-to-noise better than 10 per resolution ele- 
ment at the central wavelengths. For the purpose of this 
study, the spectroscopic data are used only for obtain- 
ing galaxy redshifts. The SDSS spectroscopic reduction 
pipeline (Schlegel et al in preparation) produces a one- 
dimensional spectrum corresponding to each fiber. The 
SDSS has an official pipeline which identifies redshifts 
(SubbaRao et al in preparation). However, the redshifts 
used in the current analysis are determined independently 
using a separate pipeline (Schlegel et al in preparation) 
whose results are nearly identical (at the 99% level) for the 
main galaxy sample to the official SDSS pipeline (which 
currently performs better for objects with unusual spectra, 
such as certain types of stars and QSOs). 



2.3. Correction to a fixed redshift frame 

Galaxies observed at different redshifts are observed 
through bandpassos at diffcrc;nt wavelengths in the 
galaxy's rest frame. The photometric measurements have 
been /C-corrected to a fixed wavelength frame using a K- 
correction pipeline (kcorrect version vl_3) described in 
detail elsewhere (Blanton et al 2002b). The method is, 
briefly, as follows. A linear combination of four SED tem- 
plates is fit to the set of broad-band magnitudes for each 
galaxy. The four templates are chosen to yield good fits 
to the broad-band magnitudes of the entire set of galax- 
ies in the SDSS spectroscopic sample. The fit provides a 
correction from the rest-frame band at the redshift of the 
galaxy to the rest-frame band at 2: = 0.1. Correction to 
the (roughly) median redshift z = 0.1 minimizes the typi- 
cal correction magnitude. The resulting fixed-frame mag- 
nitude system is denoted "° \t, °-^g, °-V, °-H, °-'^z" . The 
effective wavelengths are roughly 3220, 4340, 5660, 6940 
and 8300 A. The zeropoints are as close to AB (Oke & 
Gunn 1983) as is currently possible with SDSS calibration; 
a hypothetical source with f^, = 3.63 x lO^^'' Wm^^ Hz^^ 
at all frequencies ly would be close to mag in all five 
bands. 



2.4. Sample selection 

The sample used here is a subsample (known as 
"samples" within the SDSS collaboration) of spectro- 
scopic targets selected from all SDSS imaging targeted 
before 2001 June 20 and spectroscopically observed be- 
fore 2001 October 20. A galaxy is included in our analysis 
for this paper if it is a main sample galaxy (Strauss et 
al 2002b) satisfying the following cuts: (1) apparent mag- 
nitude 14.5 < rpet < 17.77 mag (the exact faint limit 
is variable at the 0.2 mag level, depending on the SDSS 
Main Sample photometric selection in that region of sky, 
which has varied somewhat as the survey has progressed); 
(2) mean surface brightness fj,r < 24.5 mag in 1 arcsec^ 
within the half-light radius (variable at the 0.5 mag level, 
depending again on SDSS Main Sample selection); (3) red- 
shift 0.02 < z < 0.22; and (4) fixed-frame absolute mag- 
nitude —23.5 < Mo V < —17.5 mag, corresponding to 
0.05 < (Lo.v/i*) < 12 (Blanton et al 2001). 

These cuts left a total sample of 77,461 galaxies. 

After apparent magnitude, the two most important 
cuts in our sample are on surface-brightness and abso- 
lute magnitude. It has been demonstrated that unless 
the shapes of the luminosity and surface-brightness func- 
tions change dramatically at low surface-brightnesses, the 
surface-brightness cuts exclude at most 5 percent of the 
luminosity density (Blanton et al 2001). The Schechter 
function fit to the SDSS luminosity function (from Blanton 
et al 2001 but corrected to the " V-band) has 90 percent 
of the luminosity density in the range —23.5 < Mo.i,. < 
— 17.5 mag. At least 85 percent of the luminosity density 
of the Universe is therefore in our sample; it can be used 
to make reliable estimates of the fraction of the cosmic 
luminosity density found in galaxies of different types. 



2.5. Density contributions 
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To compiitc number and luminosity densities, it is neces- 
sary to compute the number-density contribution 1/Knax 
for each galaxy, where V^ax is the volume covered by the 
survey in which this galaxy could have been observed, ac- 
counting for the flux, surface brightness, and rcdshift lim- 
its as a function of angle (Schmidt 1968). This volume is 
calculated as follows: 

1 f /•^mi„(e,« fjy 

v_=^ydo/(M)y d.-, (1) 

where f{0,4>) is described below and Zmin(^, </>) and 
Zmaxi^, 4') are defined by: 

(j)) = max(z„i^niin(^, 0), •2:^,min(^, 0), 0.02) 

{9, cf)) = min(z„i {9, (p), 0.22) . (2) 

The flux limits TOr,min(^', 4>) (equal to 14.5 mag across 
the survey) and mr,max(^, </>) (approximately 17.77 mag, 
but slightly variable as noted above) implicitly set 

m,,{„i„}(0, (/.) = Mo.v + DM(2:„,|„.„|(^, (/.)) 

+i^o.v(^„,{g;in}(^,.^)). (3) 

The surface brightness limits /ir.min(^, </') (in practice, too 
high surface-brightness to matter for this sample) and 
Mr,max(^, ^) (approximately 24.5 mag but variable in a 
known way across the survey) implicitly set Zi^^^miny{9,(j)) 
by: 

'"^.{Skx}(^' </>) = /^r- + 101ogio(l + Z^,{m„}(6',</>)) 

+Ko.^r{z^^{^^.y{e,<i>)) . (4) 

In practice, the surface brightness limits only rarely affect 
the Vinax determination. 

The function f{9,(p) is the SDSS sampling fraction of 
galaxies as a function of position on the sky. The total 
sampling rate of galaxies is computed separately for each 
region covered by a unique set of survey tiles. The nomen- 
clature of the 2dF (Percival et al 2001) is adopted; each 
such region is a "sector" (which corresponds identically an 
"overlap region" in Blanton et al 2001). That is, in the 
case of two overlapping tiles, the sampling is calculated 
separately in three sectors: the sector covered only by the 
first tile, the sector covered only by the second tile, and 
the sector covered by both. Each position on the sky is 
thereby assigned a sampling rate f{9, (j)) equal to the num- 
ber of redshifts of galaxy targets obtained in the region 
divided by the number of galaxy targets in the region. In 
regions covered by a single tile, typically 0.85 < / < 0.9; 
in multiple plate regions, typically / > 0.95. These com- 
pletenesses average to 0.92. 

3. DISTRIBUTION IN "PHOTOMETRIC SPACE" 

In this Section, the distribution of galaxies in a measure- 
ment space made up of colors, surface brightness, and con- 
centration of the radial profile is displayed and discussed. 
For the purposes of this section, the surface brightness is 
defined to be (—2.5 log^Q of) half of the Petrosian flux in 
the ° V band, divided by tt^^Iq, and the concentration c is 
defined to be the ratio ^go/^so- 

3.1. Number densities 

For any particular subset of galaxies chosen to have spe- 
cific physical properties, the sum of the inverse selection 



volumes 1/Vmax for the members of the subset is a mea- 
sure of the number density of galaxies with those physi- 
cal properties. Figure 1 shows projections of the number 
density as a function of position in the four-dimensional 
photometric space spanned by fixed-frame color "^'^(^i — g), 
fixed- frame color ^'\g — r), fixed- frame surface brightness 
fjbo.i^ and concentration 1/c. 

There are two separated density maxima apparent in 
Figure 1, one corresponding to red, concentrated, high 
surface-brightness galaxies, and one corresponding to 
bluer, less concentrated, lower surface-brightness galaxies. 
For clarity, a four-dimensional box containing the redder 
peak has been superimposed on Figure 1. The box was 
chosen by hand to be restrictive but still encompass most 
of the red peak. The box is defined by the limits: 

1.40 < ° < 2.50 

0.75 < ° H5-r) < 1.20 

18.7 < A^o-v < 20.7 

0.27 < (1/c) < 0.40 (5) 

where the colors are measured in mag and the central 
surface brightness /io.v is measured in units of mag in 
one arcsec^. We refer to the galaxies inside this four- 
dimensional box, ie, galaxies that simultaneously satisfy 
all four photometric CTits, as "red-type" galaxies. The 
separation of red-type and non-red-type galaxies is very 
similar to that shown previously to be highly correlated 
with galaxy morphological type (Strateva et al 2001). 

The integrated number density (ie, total of 1/V^ax 
values) of red-type galaxies (ie, galaxies in the four- 
dimensional box) in the sample is 20 percent of the total 
number density in samples. 

3.2. Luminosity densities 

Every galaxy in the sample has a density contribu- 
tion 1/T4iax and an °-^-band luminosity Lo.xi, computed 
with °'^pet and the luminosity distance (eg, Hogg 1999), 
so every galaxy has a luminosity density contribution 
.^^o.ij/14iax- Figure 2 shows the total luminosity density 
of the sample as a function of position in the photometric 
space. 

The integrated luminosity density of red-type galaxies 
(ie, galaxies in the four- dimensional box) in the sample is 
38 percent of the total luminosity density. 

4. DISCUSSION 

We have identified a class of "red-type" galaxies in the 
Sloan Digital Sky Survey on the basis of their colors, sur- 
face brightnesses, and radial profiles. We have shown that 
roughly two fifths of the total cosmic galaxy luminosity 
density in the visual is produced by these red galaxies. 

The high central surface brigtnesses and strong radial 
concentrations of the red galaxies suggest that they are 
red by dint of age rather than dust extinction. Galaxies 
outside the red-type box will not enter the box by becom- 
ing dustier, because the dust that makes them redder will 
also further reduce their surface brightnesses. 

It is conventional (though unnecessary in the present 
work) to type galaxies morphologically. A sample of 100 
galaxies at z < 0.1 was chosen from the red- type sample 
for morphological classification. The galaxies were chosen 
randomly, with probability proportional to number den- 
sity contribution l/Vmax (re-computed in the face of the 
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additional z < 0.1 cut). Throe authors (MB, DWH, IS) 
independently classified these 100 galaxies into the coarse 
morphological sequence E, SO, Sa, Sb/Sc, Sd/Irr. Clas- 
sification was based on ly /,y-wcighted averages of the r 
and i-band images, viewed in black and white with vari- 
able stretch on computer screens. For 95 percent of the 
sample, each pair of classifiers agree to within one mor- 
phological class in the crude sample. Across all classifiers, 
65 to 77 percent of the number-density weighted sample is 
in the E and SO classes. At most 5 percent of the sample 
is later than Sa in morphological type. Note that none 
of our conclusions depend on this by-eye morphological 
classification. 

That red galaxies contribute a large fraction of the lumi- 
nosity density is remarkable because bluer stellar popula- 
tions produce much more light per unit mass than redder. 
Typical red- type galaxies in the sample have ^-^u—g) = 
1.8 mag and ''* \g — r) = 0.9 mag; while typical galaxies 
outside the red-type region have '^■\u — g) = 1.0 mag and 
'^ \g — r) — 0.5 mag (recall that these are not zero-redshift 
colors). Conventional stellar population models suggest 
that these two colors differ in stellar mass-to- light ratio by 
a factor of two or more (Bruzual & Chariot 1993; Kauff- 
mann et al 2002), with the specific value depending on the 
assumed star formation history in each case. (Patchy or in- 
homogeneous dust extinction can modify the specific mass- 
to-light ratios, but dust obeys the same general trend, ie, 
that redder populations have higher mass-to-light ratios 
than blue.) The natural interpretation is that more than 
50 percent of the stars in the Universe are in red galaxies. 

This result is conservative, because we have only consid- 
ered a very restricted set of red-type galaxies: those that 
simultaneously satisfy all four of our cuts on color, surface 
brightness and concentration. There are certainly signif- 
icant red stellar populations contained in the centers of 
otherwise blue spiral galaxies, and old stellar populations 
contribute significantly to the stellar masses of almost all 
galaxies in the z < 0.25 Universe. A very large fraction of 
all stars are in old stellar populations. 

In conventional models of stellar population evolution, 
'^ \u — g) = 1.8 mag and roughly Solar metallicity corre- 
sponds to an age greater than about 7 Gyr (eg, Bruzual 
& Chariot 1993), or a formation redshift of 2;form > 1 (if 
the Universe is 14 Gyr old and the median redshift of the 
sample is 2; = 0.1). Such derived ages depend strongly on 



metallicity, and red galaxies are known to have abundance 
ratios different from Solar (they are alpha-enhanced; eg, 
Worthey et al 1992), so it is dangerous to take absolute 
measures of stellar population ages seriously. However, 
taken at face value, it suggests that the cosmic star for- 
mation rate was extremely high at early times. Although 
much is known about star formation history since z ^ 1 
(see references in Hogg 2001), little is known about z > 1. 
In a minimal model, the comoving volume-averaged star 
formation rate at z < 1 can be taken to go as (1 + z)'^'^ 
(Hogg 2001) and at z > 1 can be imagined to be roughly 
constant and equal to the rate at 2; ~ 1, then roughly 
65 percent of stars would be expected to be older than 
8 Gyr today. This minimal model is consistent with our 
results. Given the large munbers of old stars residing in 
red galaxies in the SDSS, supplemented by the old stars re- 
siding in bluer galaxies, it is hard see how z > 1 comoving 
star formation rate densities could have been significantly 
below 2; ~ 1 values for any significant fraction of cosmic 
time. 
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Fig. 1. — The distribution of galaxy number density in the four-dimensional photometric space of fixed-frame color ^■\u—g), fixed-frame 
color ^ \g — r), central surface brightness ixo.i^, and inverse concentration 1/c. All four photometric properties and their measurement are 
described in the text. Recall that the fixed- frame colors are fixed at redshift 0.1 not 0. The four panels on the diagonal show the one- 
dimensional distributions of each of the four photometric properties. Each histogram bin shows the total number density of galaxies in 
that bin, computed with the selection volumes Vinax- The twelve off-diagonal panels show the two-dimensional distributions for each pair of 
parameters. The shading of each pixel is proportional to the number density of galaxies with properties in that pixel, computed with the 
selection volumes Knax- Projections of the four-dimensional rectangular box that encloses the red galaxy population are superimposed on the 
off-diagonal panels. 
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